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Abstract 



We develop a non-equilibrium many-body theory of the coherent femtosec- 
ond nonlinear optical response of the Fermi edge singularity. We study the role 
of the dynamical Fermi sea response in the time-evolution of the pump-probe 
spectra. The electron-hole correlations are treated nonperturbatively with the 
time-dependent coupled cluster expansion combined with the effective Hamil- 
tonian approach. For short pulse durations, we find a non-exponential decay 
of the differential transmission during negative time delays, which is governed 
by the interactions. This is in contrast to the results obtained within the 
Hartree-Fock approximation, which predicts an exponential decay governed 
by the dephasing time. We discuss the role of the optically-induced dephasing 
effects in the coherent regime. 

PACS numbers: 71.10.Ca, 71.45.-d, 78.20.Bh, 78.47.+p 
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I. INTRODUCTION 



Ultrafast nonlinear spectroscopy provides unique and powerful tools for studying the 
dynamics of Coulomb correlation effects in semiconductors, since the time resolution can 
now be much smaller than the scattering times of the elementary excitations and, e.g, the 
period of optical phonons or plasmons. During the past decade, femtosecond pump-probe 
and wave-mixing experiments have demonstrated that in these materials both coherent 
and dissipative processes are governed by many-body effects.GJ Parallel theoretical advances 
have accounted for the two-particle correlation effects in undoped semiconductors through 
the time dependent Hartree-Fock approximation (HFA), which has developed into the very 
successful Semiconductor Bloch Equations (SBE's) formalism.!^ However, more recent ex- 
perimental observations could not be explained within this mean field approach and were 
attributed to four-particle and higher order Coulomb correlation effects. ETD These high order 
correlation effects have been mostly studied with equation-of-motion methods for density 
matricesUll or pair-operatorslll and by using the Keldysh Green's function formalism.lll 

At the same time, much less is known about the role of coherent many-body effects in the 
nonlinear optical response of doped semiconductors and, in particular, modulation-doped 
heterostructures. Some aspects of the effects of the electron-electron (e-e) scattering have 
been described within the dephasing and relaxation time approximations without micro- 
scopic justification. Most importantly however, the role of the time-dependent electron- 
hole (e-/i) correlations between the photoexcited holes and the Fermi sea (FS) electrons on 
the ultrafast dynamics is still poorly understood. The main difficulty comes from the non- 
perturbative nature of the Fermi edge singularity (FES) which dominates the absorption 
spectrum close to the absorption onset, and which cannot be described within the above 
approximat ions . Important here is that the e-h interactions between a heavy photoex- 
cited hole and the FS-electrons lead to the scattering of a macroscopic number of low-energy 
FS-pairs, which readjusts the entire FS into a new orthogonal scattering state in the course 
of the optical excitation (Anderson orthogonality catastropheE^I). 

The FES is a many-body resonance that has been observed both in doped semiconductors 
and in metals.00lii Its non-Lorentzian lineshape can be viewed as originating from the 
decay of an excitonic bound state caused by its interactions with the gapless FS excitations.0 
In a first approximation, the excitonic effects in a doped semiconductor can be described 
by extending the mean field apDroach (HFA) to include the effects of the screening and the 
Pauli blocking due to the FS.Iljii3 However, such a static treatment of the FS leads to a 
spurious bound state with respect to the Fermi energy, Ep, referred to in the following as 
the HFA bound state.lli This discrete state, with binding energy Em, would appear at the 
energy Ep — Em- Obviously, for Em < Ep (as in typical modulation-doped quantum wells), 
such a state cannot exist since it overlaps with the FS continuum, with which it interacts via 
the e-h potential.0lli The "unbinding" of this HFA bound state occurs via its interactions 
with the FS excitations,lli which are not taken into account in the HFA. Note that, in 
two-dimensional systems, a static FS allows for bound states even for arbitrarily small 
attractive interactions and therefore such an unbinding cannot arise from static screening 
or from Pauli blocking effects. This spurious bound state could be artificially merged with 
the continuum by introducing a dephasing time comparable to its binding energy. Such 
an approximation, however, neglects competely the interactions between the photoexcited 
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e-h pair and the FS excitations. A microscopic description of the unbinding of the HFA 
bound state presents a nontrivial problem due to the non-perturbative nature of the e- 
h correlations between the photoexcited hole and the FS excitations.lll Even within Green 
function techniques, this problem is rather involved because vertex correction diagrams with 
arbitrarily many crossed e-h interaction lines are as divergent as the ladder diagrams and 
should be treated on equal footing.0 To perform such a task, one must sum up at least the 
parquet diagrams and address the three-body correlations between the photoexcited hole 
and a FS excitation. Therefore, alternative methods were developed for the case of linear 
absorption, based on Fermi's golden rule with many-electron eigenstates expressed in terms 
of Slater determinants.S Such approaches become exact in the limit nie/mh ^ 1 (^e and 
nih are the electron and hole masses respectively) and describe quite accurately the FES 
lineshapes observed in typical modulation-doped quantum wells.Bil^ 

Another approach to the FES problem is based on the coupled cluster expansion 
(CCE).iiii This general many-body technique provided exact results in the limit mg/m^ -C 
10 and was used to treat the hole recoil correlations in one dimension.!^ In the latter case, 
an exact solution was obtained for uie = nth. The CCE has also been used to describe 
the e-e correlation effectH. More importantly, however, this method is well-suited for 
describing correlations in non-equilibrium systems, where it retains the advantages of dia- 
grammatic expansions without resorting to the Kadanoff-Baym ansatz or to the Markovian 
approximat ion .0 

Let us now briefly summarize recent experimental studies of dephasing and relaxation 
processes in modulation-doped quantum wells. Wang et a/.0 performed pump-probe mea- 
surements under resonant excitation conditions and observed an incoherent redshift of the 
FES at positive time delays due to the incoherent bandgap renormalization. At room tem- 
perature, Knox et. a/.S measured very fast e-e scattering times estimated ~ 10 fs, which 
implies extremely fast dephasing and relaxation. At low temperatures, Kim et a/.il measured 
long dephasing times within the frequencv range uj < Ep (the frequency domain of the FES), 
consistent with measurements in metalsEj. This result is consistent with the Landau Fermi 
liquid theoryii (which applies both to metals and doped-quantum wells due to the similar 
values of r^) and points out that, unlike for the e-h correlations, the e-e scattering is sup- 
pressed within the frequency range of the FES. Furthermore, the above experiment suggests 
that the hole dephasing times are rather long. The above experimental results were inter- 
preted within a two-level system approximation (which neglects all FES correlation effects) 
by introducing energy-dependent dephasing times obtained from the quasi-equilibrium e- 
e self-energies.0 Note here that, unlike in undoped semiconductors, the relaxation of real 
photo-excited e-h pairs due to their interactions with the FS electrons!^ can obscure the 
role of the e-h correlations. However, these incoherent effects can be suppressed by tuning 
the pump frequency below the FES resonance and thus exciting only virtual carriers. In 
this case, coherent effects will dominate. Brener et. tt/.0'il performed such pump-probe 
measurements on doped QWs with well-below-resonance pump excitation and observed 
qualitatively different nonlinear optical dynamics for the FES as compared to excitons in 
undoped QWs. This suggests that the differences in the nature of the two resonances do 
manifest themselves in the transient changes of the optical properties. 

In this paper, we develop a microscopic many-body theory for treating the coherent 
nonlinear optical response of the FES. In particular, we focus on the manifestations of the 
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non-perturbative e-h correlations (dynamical FS response) on the ultrafast time evolution 
of the pump-probe spectrum. In order to describe the dynamical FS response, we use a gen- 
eral formalism developed recently!! (which we summarize in Appendix A). Within this ap- 
proach, the transient pump-probe nonlinear absorption spectrum is mapped onto the linear 
absorption of a "pump-dressed" system, described by a time-dependent effective Hamilto- 
nian. This allows us to treat the non-equilibrium e-h correlation effects non-perturbatively 
with the time-dependent version of the CCE. We mostly focus on negative time delays and 
off-resonant excitation conditions, in which case the coherent effects dominate. 

In doped semiconductors, the pair-pair interactions, analogous to the exciton-exciton 
interactions in undoped semiconductors, are screened out, while the e-e scattering is sup- 
pressed for below-resonant excitation as dictated by Fermi Liquid theory. Instead, it is 
the e-h interactions between the photoexcited e-h pairs and the FS carriers that strongly 
affect the coherent ultrafast dynamics of the FES. An important difference from undoped 
semiconductors is that, due to its gapless pair excitation spectrum, a FS responds unadiabat- 
ically to time-dependent perturbations. In contrast, because of its finite Coulomb binding 
energy, an exciton can be polarized by the pump optical field without being ionized. A 
non-equilibrium treatment beyond the HFA is necessary in order to take into account the 
unadiabatic time-dependent change in the e-h pair-FS interactions and the e-h scattering 
processes induced by the ultrafast pump excitation. We study in detail the qualitative dif- 
ferences in the ultrafast dynamics with and without (HFA) the dynamical FS response and 
point out the important role of the e-h correlations in the coherent regime. 

The outline of this paper is as follows. In Section II we summarize our formalism for 
the nonlinear optical response and in Section III we outline the basics of the CCE approach 
to correlations. In Section IV we obtain the time-dependent parameters of the effective 
semiconductor Hamiltonian and the effective optical transition matrix elements. In Section 
V we describe the effects of the optical excitation and the many-body correlations on the 
dynamics of the photoexcited e-h pair that determines the FES lineshape. In Section VI 
we present our results for the nonlinear absorption of the FES and its time-evolution, and 
discuss the role of the dynamical FS response. Section VII concludes the paper. 

II. TRANSIENT NONLINEAR ABSORPTION SPECTRUM 

We start bvrecapping the main points of our formalism^l (see Appendix A). In the 
rotating frame,E9 the Hamiltonian describing this system is 

Htotit) = H + Hp{t) + Hs{t) . (1) 

The first term is the "bare" Hamiltonian, 

^ = E <^-q&-q + E« + ^)<«q + Vee + V^h + K,, (2) 

q q 

where aj^ is the creation operator of a conduction electron with energy and mass me, 

6Lq is the creation operator of a valence hole with energy and mass nih, Ke, K/i, and Vhh 
describe the e-e, e-h, and h-h interactions, respectively, and VL = Eg + Ep{l + me/mh)—ujp is 
the detuning of the central frequency of the optical fields, Up, from the Fermi level. Eg being 
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the bandgap (we set h = 1). The second and third terms describe the couphng of the pump 
optical field, Sp(t)e'^p'''~'''^p\ and the probe optical field, £s(t)e'^''''~"^p^*~'^\ respectively: 

Hp{t) = fiSp{t) [e'^^-^'U^ + H.c] , 

where the pump amplitude Sp{t) is centered at time t = and the probe amplitude Ss{t) is 
centered at the time delay t = t, fi is the interband transition matrix element, and 

U' = T.<bl^ (4) 
q 

is the optical transition operator. In Appendix B we clarify our convention for the time 
delay r and relate it to the most commonly used conventions in pump/probe and four wave 
mixing (FWM) experiments. 

In many experiments, the amplitude of the probe field is much smaller than that of 
the pump, ^ |^^s(^)|- In that case, as was shown in Ref. |39| (see Appendix A), the 

experimentally-measurable nonlinear optical polarization can be obtained in terms of the lin- 
ear response functions of a "pump-dressed" semiconductor to a probe field (within x^^^ , this 
is true even for comparable pulse amplitudes). This "dressed" system is described by a time- 
dependent effective Hamiltonian H{t), which is obtained by performing a time-dependent 
Schrieffer- Wolff /Van VleckiS canonical transformation on the Hamiltonian H + Hp{t). As 
we shall see in the next Section, in all the cases of interest, the effective Hamiltonian H{t) has 
the same operator form as the bare Hamiltonian H, with the important difference that the 
band dispersions (effective masses) and interaction potentials are time-dependent through 
Sp{t). Thus, the calculation of the nonlinear absorption spectrum reduces to that of the 
linear absorption spectrum of the "pump-dressed" semiconductor with uncoupled "effective 
bands" — a great simplification that allows us to use straightforward generalizations of well 
established theoretical tools. It is important to note that, as was shown in Ref. ^ (see also 
Appendix A), such a "dressed semiconductor" approach is not restricted to monochromatic 
pulses and is valid for any pulse duration. In Appendix A, we show that the pump-probe 
nonlinear polarization has the following form : 

a.(t) = -^;,2g^k..r-..,(*-.) f' dt'S,{t'){%{t)\U{t)}C{t,t')U^{t')\%{t')) . (5) 

J —oo 

Here, |$o(^)) is the state evolved with H{t) from the semiconductor ground state |0) of H, 
W{t) is the effective optical transition operator describing the probability amplitude for the 
photoexcitation of an e-h pair by the probe field in the presence of the pump excitation, 
and }C{t,t') is the time-evolution operator satisfying the Schrodinger equation 

z^/C(t,0 = i^(t)/C(t,0- (6) 

The above equation describes the time evolution of a probe-photoexcited e-h pair in the 
presence of the pump excitation. The effective Hamiltonian and effective transition operator 
are given by (see Appendix A) 
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H{t) = Ho-^ {£,it) [P(t),f/t 



H.c, 



(7) 



and 



u\t) = u^ + ^[v{t),[u\v\t) 



v\t)Au\v{t) 



(8) 



where the operator 'P^(t), which generates the canonical transformation, satisfies the equa- 
tion 



.dV\t) 
' dt 



H,V\t)] -ii8,{t)U\ 



(9) 



with the initial condition = before the pump arrives. Eqs. (^) and (§) include all the 
pump-induced corrections to U^if) and H{t) up to the second order in the pump optical field 
and are valid when [dp/Vtf' ^ 1 (for off-resonant excitation) or {dptpf' ^ 1 (for resonant 
excitation), where tp is the pump duration. 

It should be emphasized that, although Eq. (^ gives the effective Hamiltonian up to 
the second order in the pump field, EpSp-, the polarization expression Eq. (|^) describes the 
effects of H in a// orders. For example, as we shall see in Section IV, the pump-induced term 
in H contains self-energy corrections to the electron/hole energies, which describe (among 
other effects) the resonance blueshift due to the ac-Stark effect. Although the magnitude 
of these self-energy corrections, calculated using Eq. (|^, is quadratic in £p, the correct 
position of the resonance can only be obtained by evaluating the pump/probe polarization 
(§) nonperturbatively (beyond x''^'*); i-^-, without resorting to the expansion of the time- 
evolution operator /C(t, t') in the pump field. Importantly, the same is true when calculating 
the effects of the self-energy corrections on the e-h correlations. As we shall see in Section 
VI, such a nonperturbative (in the pump field) treatment of the non-linear response of the 
FES is crucial for the adequate description of the pump/probe spectrum at negative time 
delays. In Section V, we will describe the corresponding procedure, which accounts for 
the FS dynamical response. In contrast, the third-order polarization, x^^\ can be simply 
obtained from Eq. (^) by expanding }C{t,t') to the first order in the pump-induced term 
in H [second term in Eq. (|^)]. We did not include in Eq. (|^) the biexcitonic contribution 
(coming from the excitation of two e-h pairs by the pump and the probe pulses) since it 
vanishes for the negative time delays (r < 0) of interest here.ii 

The important advantage of Eq. (^, as compared to the equations of motion for the 
polarization, comes from its similarity to the linear polarization that determines the linear 
absorption spectrumll. This can be seen by setting Sp{t) = in the Eqs. (|^) and (||), in 
which case the effective time-evolution and optical transition operators transform into their 
"bare" counterparts: IC{t,t') — *■ e~*-'^'^*~* ^ and f/^(t) . Moreover, like U\ the effective 

transition operator f/^(t) creates a single e-h pair, while, as we shall see below, the effective 
Hamiltonian H{t) can be cast in a form similar to H . This allows one to interpret the 
Fourier transform of Eq. (^ as the linear absorption spectrum of a "pump-dressed" semi- 
conductor with two uncoupled but time-dependent effective bands. This mapping simplifies 
significantly the calculation of the FES ultrafast nonlinear optical response by allowing a 
straightforward generalization of the CCE. 
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III. COUPLED CLUSTER EXPANSION 



In this Section, we show how the time-dependent CCE can be used to study the ef- 
fects of the e-h correlations (dynamical FS response) on the time evolution of the e-h 
pair photoexcited by the probe. Our goal is to evaluate the many-body state = 
]C(t,t')W (t')\^Q(t')) that enters in Eq. (^. This state satisfies the Schrodinger equation 

t^^\^{t)) = m\m)- (10) 

As already mentioned, H{t) has the same form as the bare Hamiltonian H. This al- 
lows us to obtain |^E'(t)) through a straightforward generalization of the linear absorption 
calculationsiUHi. After eliminating the valence hole degrees of freedom,@'ii {"^{t)) is ex- 
pressed in the CCE form 

|^(t)) = e^W|<l>(t)), (11) 
where the time-dependent operator S{t) creates FS-pairs and is given by 



S{t) = s(p,k,t)aj,ak + ^ S2(Pi, P2, ki, ka, t)aj,^aj,^ak2aki + 

p>kp,k<kp pi,p2>kp,ki,k2<kp 



(12) 



while the state |$(t)), discussed in Section V, describes the time evolution of the probe- 
induced e-h pair. In Eq. (^2|), the amplitude s(p, k, t) describes the e-h correlations which, in 
particular, are responsible for the unbinding of the HFA bound state; the two-pair amplitude 
S2 describes the e-e interaction effects at the RPA level and beyond. iFiom a physical point 
of view, the operator e"^*^*) describes the readjustment of the FS density profile during the 
optical transition in response to the FS interactions with the photoexcited e-h pair. 

Substituting Eq. ([TT| ) into the Schrodinger equation Eq. (p!0|) , multiplying by the operator 
g-5(i) f^Qj^ ^Yie Ihs, and using the fact that [S{t), S{t')] = [S{t), S{t')] = 0, we obtain 

^^^\^t)) + ^S{tmt)) = e-^W#(t)e^W|$(t)), (13) 

where the transformed Hamiltonian on the rhs can be expressed in terms of the commutator 
series (Baker-Campbell-Hausdorff expansion) 

e-^W^(t)e^W = H{t) + [Hit), S{t)] + \[[H, S{t)i S{t)] + ■■■ (14) 
An important advantage of the CCE is that, due to the FS momentum restrictions in Eq. 



(p!^), the above series terminates after the first few terms (three for quartic interactions) and a 
closed-form expression of the transformed Hamiltonian (0) can be obtained in terms of S{t). 
By requiring that all FS-pair creation processes are eliminated from the above equation, we 
obtain the CCE equationsSS for S{t). Before proceeding with such a calculation however, 
we derive in the next Section explicit expressions for H{t) and U{t). 
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IV. EFFECTIVE HAMILTONIAN AND OPTICAL TRANSITION OPERATOR 



In this Section, we derive second-quantized expressions for the effective Hamiltonian 
H(t) and optical transition operator U{t). We start with Eqs. (|^ and (||), which express 
H{t) and U{t) in terms of the canonical transformation operator V\t). The latter is given 
by Eq. (j^), which includes the effects of the Coulomb interactions on the pump photoexcita- 
tion. It is important to realize that the effect of Coulomb e-h interactions on the pump and 
the probe photoexcitations is very different. For an adequate description of the FES, the 
e-h interactions should be taken into account non-perturbatively for the profee-photoexcited 
pair. Indeed, the nonlinear absorption spectrum at a given frequency uj close to the FES 
resonance is determined by the time-evolution for long times (of the order of the dephasing 
time T2) of an e-h pair photoexcited by the probe at the energy Since the charac- 
teristic "e-/i interaction time" i?^/ (inverse HFA bound state energy) that determines the 
nonexponential polarization decay is much shorter that the dephasing time, the long-time 
asymptotics of the response function (to the probe) depends non-perturbatively on the e-h 
interactions. In contrast, a short pump optical pulse excites a wavepacket of continuum e-h 
pair states (unlike in the discrete exciton case) with energy width ~ t~^, which thus evolves 
during timescales comparable to the pulse duration tp. Also, the corrections to the effective 
Hamiltonian are determined by the time evolution of the pump-induced carriers only up 
to times ~ tp (see Eq. (|^). Therefore, if the "e-/i interaction time" E]^ is larger than 
tp, tpEM ~ 1, (i.e., if the pump pulse frequency width exceeds Em), the e-h interactions 
can be treated perturbatively when describing the time-evolution of the j^nm^^-photoexcited 
pairs. This can also be shown explicitly for the third-order nonlinear polarization. In the 
general expression for x^^\ contributions that depend on the pump are integrated over 
the width of the pump pulse; therefore, any resonant enhancement of x^^^ fhat depends on 
the pump frequency will be broadened out for sufficiently short pulses with frequency width 
that exceeds Em- In other words, when deriving the pump-renormalized parameters, one 
can treat Coulomb interactions perturbatively if the above condition is fulfilled. In fact, 
the above situation is somewhat similar to the calculation of the linear absorption spectrum 
close to the indirect transition threshold, where perturbation theory can be used.E^'iil Thus 
the above consideration applies even for long pulse durations provided that the detuning O 
exceeds Em- This is certainly the case for the off-resonant excitation conditions considered 
in Section VI. However, in order to obtain the full absorption spectrum, the time-evolution 
of the ])rofee-photoexcited pair with such effective Hamiltonian (with perturbatively calcu- 
lated time-dependent parameters) should be treated non-perturbatively, as described in the 
following Section. 

We now proceed with the derivation of the effective Hamiltonian. To lowest order in the 
interactions, 'P^(t) can be presented as 



'PHt) = J2'Pekici,t)a\bl^ + l J2 ^:/.(PP';k;t)6Lp~p,ata;,ak 
q p,p',k 

+1 E n'k(PP;k;t)a;+p,_k6lp6lp,&-k, (15) 



P,p',k 



where Veh is the probability amplitude for excitation of an e-h pair with zero momentum 
satisfying 
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dt 



where F = T2 ^ describes the dephasing due to processes not included in the Hamiltonian H 
(e.g., due to phonons). In Eq. (p!5|), 



n\(pp';k;^) 



dt'e 



, -i{t-t')i 
e V 



y? Y* p+p'— k 



X 



{t;(p - k) [P,,(p', t') - P,,(p + p' - k, t')] - (P P')} 
describes the scattering of the photoexcited e-h pair with an electron, and 

P,\(pp'; k; t) = -z / rft'e V p' " " 1^ 

J —00 

Mp - k) [Pe.(p',t') - Pe.(p + P' - k,t')] - (P P')} 



(17) 



(18) 



describes the scattering of the photoexcited e-h pair with a hole. The above expressions 
describe in the lowest order in the screened interactionil t;(p — k) the coherent pump- 
induced processes, the effects of the Hartree-Fock pair-pair and pair-FS interactions, and 
the dynamical FS response to the pump photoexcitation. 

By substituting Eq. ([T5|) into Eq. (^, we obtain the following expression for the effective 
optical transition operator: 

UHmoit)) = E Mp(t)at 6lp|0) + ] Yl Mpp,k(t)at.a;,&Lp-p'«k|0), (19) 



p>kf 



p,p'>kp,k<kp 



where the effective matrix element Mp{t) includes corrections due to phase space filling and 
Hartree-Fock interactions, and Mpp/k(t) is the probability amplitude for indirect optical 
transitions0 induced by the pump optical field, which contribute to the pump-probe po- 
larization in the second order in the interactions. The explicit expressions for Mp{t) and 
Mpp/k(t) are given in Appendix C. 

We turn now to the effective Hamiltonian H{t). After substituting Eq. ( p!5[ ) into Eq. (|^ 
we obtain that 



(20) 



where 



Pe/^(q,t)-E^eh(qq';q;^) 



is the effective conduction electron energy; 
is the effective valence hole energy; 



n/.(q,t)-E^e\(qq';q;^) 



(21) 



(22) 
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Veh{t) = Veh{<i] kk'; t)4+q«k&Lk'-q^-k', (23) 

kk'q 

is the effective e-h interaction; and 

Vee{t) = ^Y. ^k'; t)al^^al,_^a^>a^, (24) 

^ kk'q 

is the effective e-e interaction. The exphcit expressions for t;eh(q; kk'; t) and t;ee(q; kk'; t) 
are given in Appendix C. As can be seen, H{t) has the same operator form as the bare 
Hamiltonian H. However, both the effective band dispersions and the effective interaction 
potentials are now dependent on time. Note here that the above pump-induced renormahza- 
tions only last for the pulse duration tp. As discussed above, they are therefore perturbative 
in the screened interactions for tpEM ^ 1 or for f2 > Em- 

Let us first discuss the effect of the pump-induced self-energy corrections to the conduc- 
tion and valence band energies, given by the last terms in Eqs. (^) and (|2^). The dispersion 
of the effective band is shown in Fig. 1. As can be seen, the pump pulse leads to a bandgap 
increase as well as a change in the momentum dependence (band dispersion) that last as 
long as the pump pulse. The magnitude of the bandgap increase is of the order of {iiSpY /VL 
(for off-resonant excitation) and {^£p)Hp (for resonant excitation) and leads to, e.g., the 
ac-Stark blueshift. As we shall see, for pulse duration shorter than the dephasing time, it 
also leads to bleaching and gain right below the onset of absorption, analogous to the case of 
excitons or two-level systems.!'! It should be emphasized that these are coherent effects that 
should not be confused with the incoherent bandgap redshift due to the e-e interactions 
among real photoexcited carriers.iil Note also that the above bandgap renormalization is 
induced by the transverse EM-field of the laser, as compared to the usual bandgap renor- 
malization due to a longitudinal EM-Held, i.e.. Coulomb screening. The change in the band 
dispersion, whose relative magnitude is of the order of (fiSp/Q)"^ (for off-resonant excitation) 
or [nSptpY (for resonant excitation), can be viewed as an increase in the effective density 
of states and, to the first approximation, in the effective mass. This is important in doped 
semiconductors because, as we shall see later, it leads to a time-dependent enhancement of 
the e-h interactions and scattering processes with the FS electrons. 

The effective Hamiltonian H{t) also includes pump-induced corrections in the effective 
interaction potentials, determined by the pair-pair and pair-FS interactions during the 
pump photoexcitation. By expanding Eqs. ( |C3| ) and ( |C^ ) for carrier energies close to the 
Fermi surface using Eqs. ( p^Tf ) and (|l5), one can show that these corrections vanish at the 
Fermi surface; for the typical FS excitation energies Ae ~ Em that contribute to the FES, 
their order of magnitude is {fiSp Ae / Q'^)'^ (for off-resonant excitation) or {fiSpAetp)"^ (for 
resonant excitation). Thus the corrections to the interaction potentials are suppressed for 
below-resonant excitation by a factor of {EmI^Y or for short pulses by a factor of {EMtp)"^, 
as compared to the self-energy corrections. Such a suppression is due to the Pauli blocking 
effect and the screening, which leads to the vanishing of the pump-induced corrections to 
the interaction potentials at the Fermi surface. Similarly, the pump-induced indirect optical 
transition matrix elements Mpp/k(t) are suppressed by the same factor as compared to the 
direct transition matrix element Mp(t) [first term in Eq. ([T9|)], while they contribute to the 
pump-probe polarization only in the second order in the screened interactions. Therefore, 
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in the doped case, the screened Coulomb interaction leads to subdominant parameter renor- 
malizations to the effective Hamiltonian H{t) and transition operator ?7^(t) for sufficiently 
short pump pulses or for off-resonant excitation. In the following Section, we study the 
e-h dynamics caused by the pump-induced self-energies and direct transition matrix ele- 
ments which lead to the strongest nonlinearities in the case of below-resonant or short pulse 
excitation. 



V. ELECTRON HOLE DYNAMICS 



In this Section, we derive the final formulae for the nonlinear pump-probe polarization 
of the interacting system by applying the CCE to the effective Hamiltonian H{t) in order 
to treat the dynamical FS response. The CCE equation Eq. (0) contains the operator 
S(t) described by a hierarchy of coupled equations for the amplitudes s, S2, ■ ■ ■, defined by 
Eq. (|12|) . In the coherent case of interest here, the e-e scattering effects are suppressed, 
and the nonlinear absorption spectrum is dominated by the e-h interactions. This allows 
us to use the dephasing time approximation for treating the probe-induced e-e scattering 
processeSii^, in which case the above hierarchy terminates after s. Importantly, the 
e-h correlations (dynamical FS response) are still treated non-perturbatively, since they 
are determined by s.&B Then all FS-pair creation processes can be eliminated explicitely 
from the rhs of Eq. (0), leading to the following nonlinear differential equation for the 
one-FS-pair scattering amplitude ,iHHi s(p,k, t): 



.^g(p,k,t) 
dt 



el{t)-el{t) s(p,k,t) = V 



1+ <vxt) 



E 

k'<kp 



[p, k', t) 



■ (25) 



The e-h scattering processes described by the above equation are sketched in Fig. 2a. Here 
V is the s-wave componentEl of the screened interactional, approximated for simplicity by 
its value at the Fermi energy.llliiii This neglects plasmon effects, which are however small 
within the frequency range of the FES.i'il Although this approximation is standard for 
the linear absorption case, its justification for the transient spectra requires more attention. 
Indeed, the characteristic time for screening buildup is of the order where VLp is the 
typical plasma frequency corresponding to the FS.c3@0S This time is however shorter 
than the typical pump duration ~ lOOfs and dephasing time (which is of the order of ps 
near the Fermi energ so the Coulomb interactions can be considered screened also for 

the nonlinear absorption case. In Eq. (P^), we neglected the hole recoil energy contribution 



to the excitation energyHl since s\^{t) — eo(^) ~ Em due to a sufficiently heavy hole mass.ll^ 



Note that, by increasing the hole effective mass, the pump-induced holeself-energy, Eq. 
reduces the hole recoil energy and thus the corresponding broadening.li3 In real samples, the 
relaxation of the momentum conservation condition due to the disorder will also suppress 
the hole recoil broadening effects. 

lYrom. Eq. (|l^) we then easily obtain the following expression for |$(t)): 



m)) 



p>kp 



la;6t|o), 



(26) 



where is the creation operator of the zero-momentum hole state and the e-h pair wave- 
function $p(t;t') satisfies the "Wannier-like" equation of motion. 
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.d%it,t') 



dt 



p'>kp 



where 



V{p,t) = V 



1- ^(P,k',t) 

k'<kp 



(28) 



is the effective e-h potential whose time- and momentum-dependence is determined by the 
response of the FS electrons to their interactions with the probe-induced e-h pair (vertex 
corrections) [sketched schematically in Fig. 2(b), responsible for the unbinding of the HFA 
bound state], 



eAit) =V 



k'Kkp 



J2 s{p',k',t) 

p'>kp 



(29) 



is the self-energy due the to the sudden appearance of the photoexcited hole, which leads 
to non-exponentiaFpolarization decay [described by ImeA(t)] due to the Anderson orthogo- 
nality catastropheli^ and a dynamical resonance redshift [described by ReeA(t)], and Fp de- 
scribes all additional dephasing processes (due to e-e interactions, hole recoil, and phonons). 
Eq. (13) should be solved with the initial condition %{t' ,t') = Mp(t'), where Mp(t') is de- 
fined by Eq. (0). 

It is worth stressing here the analogy between Eqs. (^61) and (pTf ) and the corresponding 
problem in undoped semiconductors. Indeed, Eq. (|26| ) is the direct analog of an exciton 
state, whereas again Eq. ( pTf ) is very similar to a Wannier equation. However, Eqs. (pSf ) 
and (^) include the effects of the interactions between the probe-photoexcited e-h pair 
and the FS-excitations, and the wavefunction $p(t,t') describes the propagation of the 
photoexcited pair "dressed" by the FS excitations. Such a "dressing" is due to the dynamical 
FS response, which leads to the dynamical screening of the effective e-/i interaction Eq. (pS]). 
The time-dependence of the latter is determined by the FS scattering amplitude s(p, k, t) 
and is affected by the pump excitation as described by Eq. (p5|). One can easily verify that 
by setting s(p, k, fL= in Eq. (|27D , we recover the results of the Hartree-Fock (ladder 
diagram, static FSli^) approximation. If one neglects the nonlinear (quadratic) term in Eq. 
(^), one recovers the three-body (Fadeev) equations.El@ Note that the coupled equations 
for $p(t,t') and s(p, k, t), obtained by neglecting the multroair excitations in Eq. (|T2|), can 
be extended to include the hole recoil-induced corrections.^ 

Using Eqs. ([H]), (|T2]), and (^), we now can express the pump-probe polarization Eq. 



(^) in terms of the e-h wavefunction $ and the effective transition matrix element Mp. 
Assuming, for simplicity, a delta-function probe pulse centered at time delay r, Ss(t) = 
SsS{t — r), we obtain a simple final expression {t > t): 



a.(t) = -^/x'^.e^'^='-*"''(*-") E Mp(t)$p(t,r) 



p>kp 



(30) 



Eq. ( pUP expresses the pump-probe polarization in terms of two physically distinct contri- 
butions. First is the effective transition matrix element Mp{t), which includes the effects of 
pair-pair and pair-FS interactions and Phase space filling effects due to the pump-induced 
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carriers present during the probe photoexcitation. Second is the wavefunction $p(t) of the 
e-h pair photoexcited by the probe, whose time dependence, determined by Eqs. ( p7|) and 
(p5|), describes the formation of the absorption resonance. Despite the formal similarities, 
there are two important differences between the doped and the undoped cases. First, in 
the doped case, the time evolution of the e-h wavefunction $p(t) is strongly affected by 
the interplay between the e-h correlations and the pump-induced transient changes in the 
bandgap and band dispersion relations. As we shall see later, this can be viewed as an 
excitation-induced dephasing. Second, unlike in the undoped casej^ the pump-induced 
corrections in the effective matrix element Mp[t) are perturbative in the screened interac- 
tions if the pump detuning or the pump frequency width exceed the Coulomb energy Em- 
In the next Section, we demonstrate the role of each of these effects both analytically (for 
CW excitation) and numerically (for short pulses). 



VI. PUMP PROBE DYNAMICS 

A. Monochromatic Excitation 

In this subsection, we analyze the FES pump-probe spectrum in the case of monochro- 
matic excitation. For monochromatic pump, the theory developed in the previous Sections 
applies for pump detunings larger than the characteristic Coulomb energy, Q > Em- Close to 
the Fermi edge, the linear absorption spectrum of the FES can be approximately described 
with the analytic expressionEl'0 

a(u;) oc Af (^^y , (31) 

where J\f is the density of states, u is the optical frequency measured from the Fermi edge, 
and f3 = 26 /ir — (S/tt)'^ is the FES exponent, where S ~ tan~^(7r5f) is the s-wave phaseshift 
of the screened e-h potential evaluated at Ep, and g = AfV is the dimensionless parameter 
characterizing the scattering strength. The monochromatic pump excitation leads to a 
resonance blueshift, originating from the shift in the effective band energies (see Fig. 1), 
and to a bleaching mainly due to the Pauli blocking which reduces the effective transition 
matrix element (analogous to the dressed atom picture^). More importantly, however, the 
pump-induced change in the band dispersion increases the density of states Af close to the 
Fermi surface and thus also increases both the e-h scattering strength g and the phaseshift 
6. This, in turn, leads to an increase in the FES exponent (3 that determines the resonance 
width and lineshape. In contrast, in the case of a bound excitonic state of dimensionality D 
and Bohr radius a^, the resonance width remains unchanged, while the oscillator strength, 
oc a^^, increases by a factor ~ (1 + D Am), where Am is the pump-induced change in the 
effective massif. Such an optically-induced enhancement of the exciton strength competes 
with the bleaching due to the Pauli blocking and exciton-exciton interactions. This results 
in an almost rigid exciton blueshift, consistent with experiment0'ii and previous theoretical 
resultsH. 

However, in the case of a FES resonance, the pump-induced change in the exponent 5 
leads to a stronger oscillator strength enhancement than for a bound exciton state. Obvi- 
ously, such an enhancement cannot be described perturbatively, i.e., with an expansion in 
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terms of the optical field, since the corresponding corrections to the absorption spectrum 
diverge logarithmically at the Fermi edge. As can be seen from Eq. (31), the effect of the 
pump on the FES can be thought of as an excitation-induced dephasing that affects the 
frequency dependence of the resonance; again, this is in contrast to the case of the exciton. 
In the time domain, this also implies a memory structure related to the response-time of the 
FS-excitations. Therefore, the qualitative differences between the nonlinear optical response 
of the FES and the exciton originate from the fact that an exciton is a discrete bound state, 
while the FES is a continuum many-body resonance. The FS responds unadiahatically to 
the pump-induced change in the density of states via an increase in the e-h scattering of 
low-energy pair excitations. Such scattering processes, which determine the response of the 
Fermi sea to the hole potential in the course of the optical excitation, are responsible for the 
unbinding and broadeningof the HFA bound state.a Therefore, the pump field changes the 
broadening and dephasing effects even for below-resonant photoexcitation. On the other 
hand, due to the finite Coulomb binding energy of the exciton, the pump optical field can 
polarize such a bound state and change its Bohr radius without ionizing it. In the next 
subsection, we study the above effects in the case of short pulse excitation. 



B. Short pulse Excitation 

We now present our results for the nonlinear absorption of the FES in the case of short 
pulse excitation. The results were obtained by solving numerically the differential equations 
(^) and (|25|) , using the Runge-Kutta method, for Gaussian pulses with duration tp = 
2.0Ep^. We only considered negative time delays and focussed on below-resonant pump 
excitation in order to suppress the incoherent effects due to the e-e scattering of real pump- 
induced carrier populations with the FS. Under such excitation conditions, the coherent 
effects in which we are interested dominate, and the Coulomb-induced corrections to the 
effective parameters, discussed in Section IV, are perturbative. Our goal is to study the 
role of the dynamical FS response {e-h correlations) on the pump-probe dynamics. For 
this reason we compare the results of our theory to those of the HFA, obtained by setting 
s(p, k, t) =0 in Eq. (|27|). As mentioned above, in the latter case, the (spurious) HFA bound 
state does not interact with the FS pair excitations, even though it can merge with the 
continuum when one introduces a very short dephasing time. 

In Fig. 3, we compare the linear absorption lineshape (in the absence of pump, Sp{t) = 0) 
of the FES to the HFA (without the dynamical FS response). We use the parameter values 
g = 0.4 and F = O.lEp, which were previously used to fit the experimental spectra in 
modulation doped quantum wells.BEl For better visibility, we shifted the curves in order 
to compare their lineshapes. The linear absorption FES lineshape is consistent with that 



obtained in Refs. |I8| and On the other hand, the HFA spectrum is characterized by 
the coexistence of the bound state and a continuum contribution due to the fact that, in 
2D, a bound state exists even for an arbitrary weak attractive potential. We note that if 
one limits oneself to linear absorption, it is possible to artificially shorten the dephasing 
time T2 = r~^, mainly determined by the hole recoil effects, by taking F ~ Em- Then 
the spurious discrete state and the continuum merge, and the discrepancy between the two 
linear absorption lineshapes decreases. This trick has been used for phenomenological fits 
of linear absorption experimental data.0 Below we show, however, that in the nonlinear 
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absorption case, the differences in tlie transient spectra are significant so tliat tlie processes 
beyond HFA can be observed experimentally. 

Let us turn to the time evolution of the pump-probe spectra. In Fig. 4 we show the 
nonlinear absorption spectra calculated by including the dynamical FS response (Fig. 4(a)) 
and within the HFA (Fig. 4(b)) at a short time delay r = —tp/2. The main features of the 
spectrum are a pump-induced resonance bleaching, blueshift, and gain right below the onset 
of absorption. For off-resonant pump, these transient effects vanish for positive time delays 
after the pump is gone, and persist for negative time delays shorter than the dephasing time 
T2 = r~^. Similar features were also obtained for different values of the pump amplitude, 
duration, and detuning. They are mainly due to the broadening induced by the transient 
renormalization of the energy band dispersion [Eqs. (^if) and (|2^) ] when its duration ~ tp 



is shorter than the dephasing time (analogous to excitons and two-level systems 

We now turn to the role of the e-h correlations. In Fig. 5 we compare the differential 
transmission spectrum calculated by including the dynamical FES response or within the 
HFA for long and short negative time delays. Note that, in pump/probe spectroscopy, the 
experimentally measured differential transmission is 

Z,5T(...)^m^ ^.fe)-7-(^.-°) . (32, 

where Ts{Sp) is the transmission coefficient in the probe direction in the presence of the 
pump field Sp. In the weak signal regime, it reproduces the changes in the probe absorption 
coefficient a{uj,T): DST{u!,t) oc —Aa{uj,T). Fig. 5(a) shows the results obtained for a 
long time delay, r = — I.5T2 = —15.0Ep^ , in which case frequency domain oscillations are 
observed. These oscillations are similar to those seen in undoped semiconductors and two- 
level systems; however, their amplitude in the FES case is reduced. On the other hand, 
as shown in Fig. 5(b), for time delays comparable to the pulse duration, r = — O.IT2 = 
-tp/2 = -1.0Ep\ the main features are a blueshift and bleaching. In this case the e-h 
correlations lead to a substantially larger width and asymmetric lineshape of the differential 
transmission spectrum. This comes from the different response of the FES to the pump- 
induced dispersion renormalizations when the e-h correlations are accounted for. This is 
more clearly seen in Fig. 6(a), where we plot the magnitude of the resonance decrease, 
evaluated at the peak frequency, as a function of r. We use, of course, the same values of 
the parameters in the two calculations and yet we find that the bleaching of the FES peak 
is substantially stronger when the dynamical FS response is included than in the HFA case. 
Note that for |r| ~ F"^ the FES resonance is actually enhanced by the pump, as can be 
seen more clearly in Fig. 7. The time dependence of the resonance bleaching is strikingly 
different in the two cases. In the HFA case, \DST{uj,t)\ evaluated at the instantaneous 
peak frequency decays over a time scale |r| ~ F^^, i.e. during the dephasing time. This is 
similar to results obtained for a two-level system with the same effective parameters. On 
the other hand, the decay of \DST{u, r)| at the peak frequency is much faster when we take 
into account the e-h correlations. Note that the above results were obtained for off-resonant 
excitation. Under resonant conditions, we find that a spectral hole is produced. In Fig. 8 
we compare the resonance blueshifts, evaluated at the peak frequency, as a function of r. 
Again, a larger blueshift is predicted when the dynamical FS response is included. This 
suggests that in the experiment of Ref. where similar blueshifts were observed in two 
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quantum well samples (one modulation doped with a FES and one undoped sample with 
a 2D-exciton) the effective parameters were larger in the latter case, due to the absence of 
screening and exciton-exciton interaction effects. 

In order to gain qualitative understanding of the role of the e-h correlations, let us for 
a moment neglect the momentum dependence of the pump-renormalization of the band 
dispersion and the phase space filling effects and consider the bleaching caused by a rigid 
semiconductor band shift AEg{t), obtained from the pump-induced self-energies, Eqs. ( PT] ) 
and (0), evaluated at the bottom of the band (note that AEg{t) lasts for the duration of 
pulse). Within this model, the pump excitation has no effect on the e-h scattering amplitude 
s(p,k, t) [see Eq. (^)]. It is thus convenient to factorize the effects of the rigid band shift 
on the e-h wavefunction ^p{t,t'): 

%it, t') = e-'It' ^i^9{t")dt" ^/^ _ (33) 

This relation is general and defines $p(t,t'), which does not depend on AEg{t). In the 
special case of a rigid shift, ^p(t,t') coincides with $p(t — t') describing the propagation of 
the probe-photoexcited e-h pair in the absence of the pump pulse. By substituting into Eq. 
(0) the long-time asymptotic expression^ l>p(t, t') = ^^{t — t') oc [i{t — t')EF]'^~^ that gives 
the linear absorption spectrum of the FES at ~ Ep, and substituting the resulting $p(t, t') 
into Eq. (|30|) , we obtain a simple analytic expression for the effect of a pump-induced rigid 
band shift on the nonlinear absorption spectrum: 

a{uj) oc Re ^^e^(-+^r)(t-r)-./>i^,(Odt'[,(^ _ r)EFf-\ (34) 

For AEg{t) = one, of course, recovers the linear FES absorption in the vicinity of the 
Fermi edge.@ For /? = 0, Eq. ( |5^ gives the absorption of the non-interacting continuum. 

The physics of the FES can be seen from Eq. (|34D. For /3 = 1, this gives a discrete 
Lorentzian peak corresponding to the HFA bound state. However, during the optical tran- 
sition, the e-h pair interacts with the FS electrons, leading to the and the readjustment of 
the FS density profile via the scattering of FS pairs. This results in the broadening of the 
discrete HFA bound state, which is governed by the time evolution of the FS. Such time 
evolution is unadiabatic due to the low-energy FS pairs, which leads to the characteristic 
power-law time dependence of the broadening factor in Eq. (P^. The interaction with the 
FS-pairs determines the exponent, < /5 < 1, of the latter, which leads to a non-Lorentzian 
lineshape in the frequency domain and a non-exponential decay in the time domain. A 
detailed discussion of the above physics and the analogy to phonon sidebands and collision 
broadening may be found in Ref. ^ Here it is essential to use the CCE in order to calculate 
the spectrum at all frequencies (and not just asymptotically close to Ep as with Eq. ([3^)) 
and, most importantly, to describe the non-equilibrium FS and e-h pair response to the 
time-dependent increase in the effective mass/density of states, not included in Eq. (p^). 

The resonance bleaching obtained from Eq. (|3^ as a function of r, is shown in Fig. 
6(b) for P = 0.6, corresponding to the value of the parameters used in Fig. 6(a) {g = 0.4), 
together with the HFA result (/9 = 1). Comparing Figs. 6(a) and (b), one can see that 
the rigid band shift approximation qualitatively accounts for the dynamics, but that there 
are strong discrepancies (see vertical scales), whose origin is discussed below. Both the 
magnitude and the time-dependence of the bleaching depends critically on the value of 
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which characterizes the interaction of the photoexcited e-h pair with the FS excitations. 
Because of such couphng, many polarization components are excited in the case of the 
continuum FES resonance, and it is their interference that governs the dynamics of the 
pump/probe signal. Such interference is also responsible for the resonance enhancement 
and differential transmission oscillations at r < shown in Figs. 5 and 7. As (3 increases 
the interference effects are suppressed because the energy width of the continuum states 
contributing to the FES narrows. In fact, this energy width is directly related to that of 
the linear absorption resonance. This is clearly seen in Fig. 9 where we show the effect of 
increasing g on the dynamics of the bleaching. It becomes more bound-state-like as, with 
increasing g, the FES resonance becomes narrower. On the other hand, in the HFA case, 
the decay rate is ~ T2 = F"^, i. e. it is independent on g, when Em becomes smaller than 
F, while for Em — F, the contribution of the continuum states produce a faster decay. 

Although the transient rigid band shift approximation, Eq. (0), explains some of the 
features of the dynamics of the bleaching, it strongly overestimates its magnitude. This 
is because Eq. ( p4D neglects the response of the many-body system to the pump-induced 
renormalization of the band's dispersion, neglected in Eq. (p^. Such a transient change in 
the dispersion, which can be viewed as an increase in the density of states/effective mass 
for the duration of the pump, is important because it results in an enhancement of the e-h 
scattering. For example, in the case of monochromatic excitation, this leads to the change 
in the exponent (3 of the broadening prefactor in the integrand of Eq. (|3^), discussed in the 
previous subsection. For the short-pulse case, it is not possible to describe analytically the 
effect of the pump on the e-h scattering processes, due to the non-equilibrium unadiabatic 
FS response. The latter can be described with the numerical solution of Eqs. ( pS] ) and (p7D, 
which is presented in Figs. 10 and 11 and discussed below. 

In order to show the role of the pump-induced renormalization of the band dispersion 
in the presense of the dynamical FS response, we plot in Fig. 10 the function 

F(^,r) = Im Yl %i^,r) (35) 

p>kp 

where $p(u;,r) is the Fourier transform of <l'p(t,r) defined by Eq. (|3^). Note that, in the 
presence of the band dispersion renormalization, the wave-function $p (which is independent 
of AEg{t)) no longer coincides with $p as in Eq. (|34D. As can be seen in Fig. 10(a), when 
the e-h correlations are taken into account, the pump-induced redistribution of oscillator 
strength between the states of the continuum that contribute to the resonance manifests 
itself as a dynamical redshift. This shift opposes the rigid band blueshift AEg{t) (when 
the latter is included). At the same time, the resonance strength is enhanced significantly. 
The latter effect originates from the interplay between the transient increase in the effective 
mass/density of states of the photoexcited e-h pair and the "dressing" of this pair with 
the FS excitations [described by the effective potential V{p,t) in Eq. (0)]. In contrast, 
such an oscillator strength enhancement is suppressed in the HFA (which neglects the e-h 
correlations), as seen in Fig. 10(b), in which case the main feature is the redshift of the 
resonance due to the pump-induced increase of the binding energy Em coming from the 
transient increase in the effective mass [see Fig. 10(b)]. 

In Fig. 11 we show the effect of the renormalization of the band dispersion on the non- 
linear absorption spectrum. The optically-induced increase in the e-h interactions enhances 
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significantly the strength of the FES and compensates part of the bleaching induced by the 
rigid band shift. A smaller enhancement is also seen in the HFA, where the pump-induced 
increase in the binding energy Em competes with the effects of the bandgap renormalization. 

VII. CONCLUSIONS 

In summary, we developed a theory for the ultrafast nonlinear optical response of the 
FES. We focussed on coherent effects, which dominate the pump-probe spectra during 
negative time delays and off-resonant excitation conditions. We demonstrated that the 
dynamical FS response leads to qualitatively different coherent dynamics of the FES as 
compared to the Hartree-Fock approximation. In particular, in the former case, the time 
evolution of the resonance bleaching is governed by the dephasing time, while in the former 
case polarization interference effects dominate. This results in faster FES dynamics (which 
depends on the strength of the screened e-h potential) as well as an apparent resonance 
enhancement during negative time delays. Such dynamical features should be observable 
in ultrafast pump/probe experiments. Using a simple model, we showed that the different 
dynamics of the FES and Hartree-Fock treatment can be attributed to the non-Lorentzian 
broadening of the HFA bound state due to its interactions with the gapless FS excitations, 
a process which is, of course, beyond the dephasing time approximation. In addition, we 
showed that the pump excitation directly affects the strength of the e-h scattering processes, 
which changes the frequency dependence of the resonance. The latter can be thought of as 
an excitation-induced dephasing effect that leads to a transient enhancement of the FES. 
Our results indicate that ultrafast spectroscopy provides a powerful tool to study the role 
of correlations in the nonlinear response of a Fermi liquid to the optical excitation during 
time scales shorter than the dephasing times. The dynamical features discussed above can 
also be used as an experimental signal to probe the crossover from FES to exciton bound 
state (exciton Mott transition) as a function of the FS density. 
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APPENDIX A: 

In this appendix we briefly outline our formalism.^ The pump-probe signal is determined 
by the polarization, 

P[t) = fxe-'''^'{^{t)\U\^it)), (Al) 
where the state |^E'(t)) satisfies the time-dependent Schrodinger equation. 
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\^{t)) 



(A2) 



with the Hamiltonians H and Hp^s{t), given by Eqs. 



The Hilbert space of the 



bare semiconductor, i.e. in the absence of optical fields, consists of disconnected subspaces 
Ci'^eh} which are labeled by the number of (interband) e-h pairs, u^h- The corresponding 
bare Hamiltonian, H, conserves the number of e-h pairs in each band separately and in 
the ('{'^e/ij-basis has a block-diagonal form. The Hamiltonians Hp s{t) couple the different 
subspaces Ci^eh} by causing interband transitions. 

In the description of pump/probe experiments, we are interested only in the polarization 
component propagating along the probe direction k^. For a weak probe, the nonlinear 
signal then arises from the linear response of the pump/bare-semiconductor coupled system, 
described by the time-dependent Hamiltonian H+Hp(t), to the probe-induced perturbation 
Hs{t). Within x^^\ the above is true even for comparable pump and probe amplitudes. 
However, since, in contrast to H, the Hamiltonian H + Hp{t) does not conserve the number 
of carriers in each band separately, the calculation of the linear response function is not 
practical. Therefore, we seek to replace H + Hp{t) by an effective Hamiltonian H{t) that 
does conserve the number of e-h pairs in each of its Hilbert subspaces (i.e., is "block- 
diagonal"). As derived in Ref. this can be accomplished in any given order in the 
pump field Sp{t) and for any pulse duration by using a time-dependent Schrieffer-Wolff 
/Van Vleck canonical transformationiillii. Here it is sufficient to "block-diagonalize" the 
Hamiltonian H + Hp{t) up to the second order in £p{t). The transformation that achieves 
this has the form e~^'^e~'^^[H + Hp{t)\e^^e^'^ , where the anti-Hermitian operators Ti{t) and 
T2{t) create/annihilate one and two e-h pairs, respectively. 

We proceed with the first step and eliminate the single-pair pump-induced transitions in 
the time-dependent Schrodinger equation of the pump/bare-semiconductor system. 



Hp{t) 



(A3) 



This is achieved by substituting \^{t)) = e^^*-'|x(t)) and acting with the operator e on 
the Ihs of Eq. 



-Tiit) 



3^^^*)|x(t))=e-^^W [if,(t)]e^^W|x(t)). 



(A4) 



The anti-Hermitian operator Ti(t) has a decomposition 



f,it)=Vit)e 



(A5) 



where 'P^(t) and P(t) create and annihilate single e-h pairs, respectively. The effective 
Hamiltonian can be found from the condition that the terms describing single-pair inter- 
band transitions cancel each other in Eq. (|A^). In Ref. it was shown that the multiple 
commutators of V{t) with its time derivatives can be eliminated from Eq. ( jX^ ) to all orders. 
By expanding Eq. (|A4|) and neglecting third or higher order terms in V, we obtain the 
following equation: 
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' dt 



(A6) 



with initial condition ^^(—00) = 0. The formal solution is 

J —00 



(AT) 



Note that, since the Hamiltonian H conserves the number of e-/i pairs and the optical transi- 
tion operator creates a single e-h pair, V^{t) also creates a single e-h pair. Furthermore, 
since both H and [/^ conserve momentum, so does P^{t). Eq. ( |A^ ) then takes the form 



\x{t)) 



'2 



H.c. 



\x{t)), 



where 



H.c. 



(A8) 



(A9) 



is the sought time-dependent effective Hamiltonian that conserves the number of e-h pairs 
and V^{t) is given by Eq. ([A^). Note that, since 'P^(t) is linear in the pump field £p, the 
pump-induced term in H(t) [second term in Eq. (|A9|) ] is quadratic (oc Sp£*). The rhs of 
Eq. ( |A8|) describes the pump-induced two-pair transitions. These can be eliminated as well 
by performing a second canonical transformation, = e^^*-*''|$(t)). Following the same 

procedure, we use the anti-Hermiticity of T2{t) to decompose it as 



f2(t) = A(t)e 



-H{t)e 



2i]<.r,r 



(AlO) 



where Vli't) and V2(t) create and annihilate two e-h pairs, respectively. Substituting the 



above expression for into Eq. ( |A8| ) and requiring that all two-pair transitions cancel 

out, we obtain the following equation for Vlit), 



dt 



(All) 



Note that Vlit) only affects the pump/probe polarization via higher order {S^) corrections, 
which are neglected here. However, it does determine the four-wave-mixing (FWM) polar- 
ization (see below). 

To obtain the condition of validity of this approach, it is useful to write down a formal 
solution ( [A7| ) of Eq. (|A6|) in the basis of the N-hole many-body eigenstates, \aN), with 
energies EaN, of the Hamiltonian H. Here the index a labels all the other quantum numbers, 
so that A^=0 corresponds to the semiconductor ground state |0), |al) denotes the one-pair 
states (exciton eigenstates in the undoped case, with a labeling both bound and scattering 
states), \a2) denotes the two-pair (biexciton in the undoped case) eigenstates, etc. In this 
basis, the solution of Eq. ( |A6| ) can be written as 



{(3N + l\V^{t)\aN) 
{(3N+l\W\aN) 



dt'Sp(t')e 



-r(t-t') 



(A12) 
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where we separated out the detuning Q and denoted AE^i^ = Ep^^i — EaN- It can be seen 
that for resonant excitation (small Q) the rhs of Eq. ( A12| ) is of the order of fiSptp. Thus, 
for short pulses, this parameter justifies the expansion in terms of the optical fields. For off- 
resonant excitation, this expansion is valid for any pulse duration provided that fiSp/Q ^ 1. 
Similar conditions can be obtained for the two-pair transition described by V2- 

The nonlinear polarization Eq. ( |Al|) can now be expressed in terms of the linear response 
to the probe field: 



P(t) = /ie-*^''*(^(t)|f/|^(t)) = iie~'^'^'mt)\UTit)\^{t)), 
where, in the first order in Ss{t), the state is given by 



= - I dt'}C{t,t') £s{ty'^^'^^^ul{t') + H.c. |$o(t'))- 



Here }C{t,t') is the time-evolution operator satisfying 



(A13) 



(A14) 



(A15) 



and U^{t) = e~'^^^^^e~^'^^^^W e^'^^^^e^'-'^^^ is the (transformed) optical transition operator. In 
Eq. ( |A14| ), |$o(^)) = ^{t, ~°o)|0) is the time-evolved ground state |0); since H(t) conserves 
the number of e-h pairs, |$o(^)) contains no interband e-h pairs (in undoped semiconduc- 
tors, it coincides with the ground state, |$o(^)) = |0)). /,From Eqs. ( |A14| ) and (|A13|) , the 
polarization P{t) takes the form 

J —oo 

{^o{t)\UTit)^t,t')[£,{t')e'^'-'''^'^''^Ul{t')+£:{t')e-'^^^^^^^ 



X 



$o(t')l ^s(t')e"'"'+'"''"t/^(t') +^s*(^')e~'''"'~'"""f/r(t') }C{t' ,t)UT{t)\%it)) 



. (A16) 



The above expression for the total polarization contains contributions propagating in various 
directions. To obtain the polarization propagating in a specific direction, one has to expand 
the effective-transition operator [/^(t) in terms of Ti and T2. Using Eqs. ( |A^ ) and ( |A10| ) 
and keeping only terms contributing to pump/probe and FWM polarizations, we obtainll 



f/t (t) = f/t(t) + f/t(t)e*k,r ^ f^^^^(^)g-2.k,r ^ . . . ^ 

where (to lowest order in the pump field) 



(A17) 



U 



f/|(t) = f/t + i[p(t),[f/t,pt(i) 

ul{t) = [v\t),u\ 

1 r 
2 



v^it),\uKvit) 



t 

FWM 



it) 



U,V\t)] ,V\t)] - \u,vl{t) 



(A18) 



Here operators Ul{t) = U{t) and Upwiviit) create one e-/i pair, while Ulit) creates two e-h 
pairs (note that Upwhi iii Eq- (|A18|) annihilates an e-/ipair). 
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Pump/probe polarization 



In order to extract the pump/probe polarization from Eq. ( |Al(j| ), one should retain only 



terms that are proportional to e*''^ '". Substituting Eqs. (|A17|) into Eq. ( |A16| ), we obtain 
a.(t) = P£^(t)+P£^(t), where 

p£)(t) = -z/x2e^k..r-^.,(*-.) r' dt'Ss{t'){Mt)mtMt,t')U\t')\%{t')), (A19) 

J — oo 



P^\t) = -^^^^^^■'^-^■Mt-r) f dt'S,{t'){%{t)\U2mMUlit')\%{t')). (A20) 



and 



Note that the above formulae apply to both undoped and doped semiconductors. 

Eqs. (|A19|) -( |A2CI| ) express the nonlinear pump/probe polarization in terms of the linear 
response to the probe field of a system described by the time-dependent effective Hamilto- 
nian (|A9|) . Such a form for the nonlinear response allows one to distinguish between two 
physically distinct contributions to the optical nonlinearities. Assuming that a short probe 
pulse arrives at t = r, consider the first term, Eq. ( |A19| ), which gives the single-pair (exciton 
for undoped case) contribution to the pump/probe polarization. At negative time delays, 
r < 0, the probe excites an e-h pair, described by state W {t)\^q{t)) ~ U^O); since the 
probe arrives before the pump, the effective transition operator coincides with the "bare" one 
[see Eqs. ( |A18| ) and (|A7|) ]. The first contribution to the optical nonlinearities comes from 



the effective Hamiltonian, H{t), governing the propagation of that interacting e-h pair in 
the interval (r, t) via the time-evolution operator /C(t,r). Note that since the pump pulse 
arrives at t = 0, for |r| ^ F"^, the negative time-delay signal vanishes. At t > 0, the 
e-h pair (exciton in the undoped case) "feels" the effect of the pump via mainly the tran- 
sient bandgap shift, leading, e.g., to ac-Stark effect, and the change in the band dispersions 
(increase in effective mass/density of states), leading to enhanced e-h scattering (exciton 
binding energy in the undoped case). Note that H{t) also contains a contribution coming 
from the interactions between probe- and pump-excited e-h pairs, which are however per- 
turbative in the doped case for short pulses or off-resonant excitation [see Section IV] and 
lead to subdominant corrections. Importantly, the response of the system to the optically- 
induced corrections in H(t) takes into account a// orders in the pump field, which is necessary 
for the adequate description, e.g., of the ac-Stark effect and the pump-induced changes in 
the e-h correlations. Indeed, although the pump-induced term in Eq. (^^) is quadratic in 
Sp, the time-evolution of the interacting e-h pair is described without expanding /C(t, r) in 
the pump field. On the other hand, the third-order polarization (x^^^) can be obtained by 
expanding ]C{t, r) to the lowest order. The second contribution to the optical nonlinearities 
comes from the matrix element of the final state, {^o{t)\U{t) in Eq. (|A19D . The latter, given 
by Eq. (A18), contains the lowest order (quadratic) pump-induced terms which describe 



the Pauli blocking, pair-pair, and pair-FS interaction effects (exciton-exciton interactions 
in the undoped casetHi). Note that the matrix element of the initial state contributes for 
positive time delays, i.e., if the probe arrives after the pump pulse. In this case, however, 
the pump-induced term in the effective Hamiltonian ( [A9| ) vanishes (since it lasts only for 
the duration of the pulse) so that for positive t > tp the pump/probe signal is determined 
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by the matrix elements rather than by H{t). If the probe arrives during the interaction of 
the system with the pump pulse (r ~ tp), both the effective Hamiltonian and the matrix 
elements contribute to the polarization. In this case, there is also a biexcitonic contribution 
[given by Eq. ( A20| )] coming from a simultaneous excitation of two e-h pairs by the pump 
and the probe. However, such a biexciton state does not contribute to negative (r < 0) time 
delays. As can be seen from the above discussion, our theory separates out a number of 
contributions that play a different role for different time delays and excitation conditions. 



FWM polarization 



By extracting from Eq. ( [A16| ) all the terms propagating in the the FWM direction. 



2kp — ks, we obtainil (for delta-function probe Ss{t) = Ss6{t — r)) 



FWM\ 



(t) = _,^2g.(2k,-k.).r-..,(i+.)^* ($o(t)|f//C(t,r)f/;^^(r)|$o(r)) - {t ^ r) 



, (A21) 



where the FWM transition operator Up^rJ^(t) is given by Eq. ( |A18| ) and t > t. It is 
convenient to express Up^yj^jit) in terms of the "irreducible" two-pair operator W'^it) = 
ipt2 _pt^ satisfying 



dWUt) 



H{t),W\t)] - fiSp{t)U^P{t). (A22) 



In terms of W^{t), the state UlwMit)\^o{t)) in Eq. ([A2l| ) can be presented as a sum of two- 
and one-pair contributions: 

UlwM{t)\Mt)) = UW\t)\Mt)) - V\t)UV\moit)). (A23) 

The operator W^{t), being quadratic in the pump field (oc E^), describes the simultaneous 
excitation of two interacting e-h pairs by the pump pulse and includes the two-photon 
coherence effects (biexciton and exciton-exciton scattering effects in the undoped case). The 
corresponding state, W'^ {t)\^Q{t)) , satisfies the four-particle (two-exciton for the undoped 
case) Schrodinger-like equation with a source term [coming from the second term in the 
rhs of Eq. ( |A22|) ] and can be expressed in terms of the two-exciton Green functionll. The 
first term in Eq. (|A23| ) is responsible for the (interaction-induced) finite FWM signal when 
the probe arrives after the pump.u The second term in Eq. ( |A23| ), after being substituted 
into Eq. (|A21|) , describes the diffraction of the pump field on the grating kp — k^ due to 
the interference of the pump and probe electric fields, and contributes to the Pauli blocking 
and single-exciton effects. Note that, similarly to the pump/probe polarization, Eq. (|A21|) 



gives contributions to the FWM polarization in all orders in the pump field; in this case, 
however, the third-order polarization x^^^ can be obtained by simply neglecting the difference 
between H{t) and H in the time-evolution operator }C{t,T) and in the equation for W^{t) 
because Up^rM{'t) is already quadratic in the pump field. For x^^\ the connection between 
our formalism and that of Axt and StahlSH can be established by considering the matrix 
elements of the operator W'^{t) between two e-h pair states.!! 
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APPENDIX B: 



In this appendix we clarify our convention for the time delay r and relate it to 
the most commonly used conventions in Pump/probe and Four Wave Mixing experi- 
ments. In the generic experimental configuration two laser pulses £'i(i)e*'*i'''~"^^*~*i\ and 
i?2(t)e*''2 '"~*'^'^*~*2\ respectively centered at time t — ti and t — t2 are incident on a sample. 
Let us define At as, 

At^ti- 12 (Bl) 

and consider a FWM experiment where the signal is measured in the direction 2k2 — ki. 
Then for a two- level-atom, the signal vanishes for At < 0, while for At > its amplitude, 
which decays with time as e'"*/-^^, is determined by the Pauli blocking. In a system with 
Coulomb interactions (such as a semiconductor) a FWM signal is observed both for At < 
and At > 0. The At < signal is entirely due to the Coulomb interaction. 

In pump/probe experiments, one usually chooses one pulse (the "pump") to have an 
amplitude Ep much larger than that of the other pulse (the "probe"). Eg. As discussed in 
the text, a weak probe measures the linear response of the system (bare or dressed by the 
pump). If we choose that Ep — E2, the pump induces coherent and incoherent populations 
when it arrives in the sample before the probe i.e. for t2 < ti. This is usually defined as 
"positive" time delay t = t2 — ti > in the pump/probe literature. Note that r = —At, 
i.e., the "regular" sequence in pump/probe experiments is the reverse of that of FWM 
experiments. For r < 0, the origin of the pump/probe signal is that the probe creates 
a linear polarization in the sample which lasts for time ~ T2 = F"^ and, consequently, 
is scattered by polarization excited by the pump field. The signal observed for r < is 
therefore due to coherent effects. 

In FWM experiments, there is no restriction on the magnitude of the two incident fields 
E2 and El, which are often chosen to have amplitudes of the same order. Note however 
that, at the x^^^ level, the FWM and pump-probe polarizations are linear in the Ei{t) field 
and thus the above linear response calculation applies even for comparable pump and probe 
amplitudes. 



APPENDIX C: 

In this Appendix we present the explicit expressions for the renormalized transition 
matrix elements in the presence of the pump excitation. The direct transition matrix element 
is given by 



Mp{t)^l-\Veh{p,t)f + 



KM E n\(pk';k';i) + H.c. 

k'<kp 



p >kp 

+ 2 ^ 

^ p'>kF L 



n\(p'p; p'; t) - n Jp'p; p'; t) - Khip'p; p; t) + n\(p'p; p; t) 



n\(pp';p;i) + n\(pp';p;i) 



(Cl) 
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The first term on the rhs of the above equation describes the phase space filhng contribution, 
while the rest of the terms are due to the mean field pair-pair and pair-FS interactions. 
The pump-induced indirect transition matrix element is given by 



Mpp,k(t) 



Veh{Kt)-Veh{p + p' -^,t) 



n\(PP;k;t) 



n\(pp'; P; t) + Pi*(pp'; p + p' - k; t) 



+Veh(p + P' -Kt) 



r:,{K P + p' - k; p'; t) - n\(k, p + p' - k; p; t) 



V^^{k, p + p' - k; p'; t) + VUp, p'; k; t) 



-KhiP: P'; P + p' - k; - V^^iK P + p' - k; p; t) 

+n/^(p, t)V!l{k, p + p' - k; p'; t) - Vehip', t)V2{k, p + p' - k; p; t). (C2) 
The effective e-h potential is given by 



Vehiq.; kk'; t) = 'u(q) - | Sp{t) 



V!^{k + q, k'; k' + q; t) + V^^{k, k' + q; k'; t) 
Vl,ik, k' + q; k'; t) + V^^ik + q, k'; k' + q; t) 



and the effective e-e potential is given by 



Vee{ci;kk';t)=v{q) + ^Sp{t) 



(C3) 



n\(k + q, k' - q; k'; t) - Vl^{k + q, k' - q; k; t) 
V:,{k, k'; k' - q; - Vt^{k, k'; k + q;t) 



(C4) 
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FIGURES 



FIG. 1. The effect of the pump pulse on the band dispersion, ek{t) = e^{t) + e'"_]^{t), of the 
"pump-dressed" system. SoUd hne: "bare" dispersion {t = — oo). Dashed line: pump-renormalized 
dispersion (t = 0). The bands become "heavier" for the duration of the pump. 

FIG. 2. (a) The e-h scattering processes that contribute to the rhs of Eq. (p5|). Full lines 
correspond to s(p,k, t) and thin lines to V. The diagrams describe (from top to bottom) (i) Born 
scattering of a FS electron, (ii) FS electron ladder diagrams, (iii) FS hole ladder diagrams, and 
(iv) Nonlinear vertex corrections due to the dynamical FS response, (b) Scattering processes state 
that determine the time- and momentum-dependence of the effective e-h potential V^(p, t) and 
lead to the unbinding of the HFA bound state. 

FIG. 3. Linear absorption resonance lineshape for FES (solid curve) compared to the HFA 
(dashed curve), calculated with g = 0.4 and F = O.lEp. The HFA resonance position was shifted 
for better visibility. 

FIG. 4. The absorption spectra for the FES (a) compared to the HFA (b). Solid curves: linear 
absorption. Dashed curves: nonlinear absorption. The curves were calculated with g = 0.4 and 
F = O.lEp at time delay r = — O.IF"^ = —tp/2 and pulse duration tp = 2.0Ep^ . The nonlinear 
absorption lineshapes exhibit bleaching, resonance blueshift, and gain below the absorption onset 
that differ in the two cases. 

FIG. 5. Differential transmission lineshape for various time delays calculated with g = 0.4, 
F = O.lEp, and tp = 2.0Ep^. (a) Long time delay, r = -LSF"! = -15.0^^^ For the FES, 
the oscillations in the differential transmission spectra are reduced (solid curve) as compared to 
the HFA (dashed curve), (b) Short time delay, r = -0.ir2 = -tp/2 = -l.OEp'^. For the FES, 
the differential transmission spectrum is asymmetric (solid curve) as compared to the symmetric 
lineshape for the HFA (dashed curve). The above curves were shifted for better visibility. 

FIG. 6. (a) Nonlinear absorption resonance bleaching evaluated at the instantaneous peak 
frequency as function of time delay for the FES (solid curve) compared to the HFA (dashed curve) . 
The curves were calculated with g = 0.4, F = O.lEp, and tp = 2.0Ep^ . The time-dependence of 
the pump pulse is also presented for comparison (dotted curve), (b) Same calculated using the 
rigid band shift model, Eq. (|34|). 

FIG. 7. The enhancement of the nonlinear absorption resonance of the FES (dashed curve) 
vs. linear absorption (solid curve) at long time delay r = — l.SF"^. The curves were calculated 
with g = 0.4, F = O.lEp, and tp = 2.QEp^ . 
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FIG. 8. Nonlinear absorption resonance blueshift as function of time delay for the FES (solid 
curve) compared to the HFA (dashed curve). The blueshift is significantly weaker for the HFA. 
The curves were calculated with g = 0.4, F = O.lEp, and tp = 2.0Ep^. The time-dependence of 
the pump pulse is also presented for comparison (dotted curve). 

FIG. 9. Resonance bleaching as function of time delay for the FES (a) compared to the HFA 
(b) for different strengths of the e-h interaction: g=0.4 (solid curve) and g=0.3 (dashed curve). 
The curves were calculated with F = O.lEp, and tp = 2.0Ep^ . The time-dependence of the pump 
pulse is also presented for comparison (dotted curve). 

FIG. 10. The effect of the pump-induced renormalization of the band dispersions on the e-h 
interactions. The function F{uj,t), given by Eq. (|35|), for the FES (a) compared to the HFA (b) 
in the presence (dashed curve) and absence (solid curve) of the pump pulse. The curves were 
calculated with g = 0.4, F = O.lEp, and tp = 2.0Ep^ . 

FIG. 11. The effect of the pump-induced renormalization of the band dispersions on the 
resonance strength. The nonlinear absorption spectrum for the FES (a) compared to the HFA (b). 
Solid curves: Linear absorption. Dashed curves: Nonlinear absorption. Dotted curves: Nonlinear 
absorption for a rigid band shift only. The curves were calculated with g = 0.4, F = O.lEp, and 
tp = 2.QEp^ . 
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FIG. 3 



33 



ABSORPTION (arb. units) 



I I 

o o o o o o 

N> N> b> N> N> b> 

I 

N> I 1 ^ ' ^ ' ^ \ 1 ^ ' ^ 



CO 




DIFF. TRANSMISSION (arb. units) 



II II 

o o o o o o o 




BLEACHING (percentage) 



CO ro -»■ I I 

aicnoioicn ocnocn 





FIG. 7 



37 



0.06 



0.05 




-6.0 -4.0 -2.0 0.0 



FIG. 8 



38 



BLEACHING (percentage) 



I I I 

-..I ^3 _L _.. 

I o cn o o o o o 



CO 
CO 




F(C0,T) 

o o o O 

o ji. bo o ji. bo k) 




ABSORPTION (arb. units) 



I 

o o o o o 

N>N>b>-^oj^C»N> 

I 

N> I I ^ ' ^ ' ^ \ ' ^ ' ^ ' 




